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We consider Laplace operators on metric graphs, networks of one-dimensional 
line segments (bonds) , with matching conditions at the vertices that make the 
operator self-adjoint. Such quantum graphs provide a simple model of quantum 
£N) ■ mechanics in a classically chaotic system with multiple scales corresponding to 

^ ' the lengths of the bonds. For graph Laplacians we briefly report results for the 

VP spectral determinant, vacuum energy and heat kernel asymptotics of general 

, graphs in terms of the vertex matching conditions. 

o 
o 



1. Introduction 

Quantum graphs provide a simple, analytically tractable, model in which 
to investigate quantum phenomena. This is of particular interest for the 
study of vacuum energy and related spectral quantities as the spectrum 
of a quantum graph generically incorporates the features associated with 
quantum mechanical systems that are classically chaotic. 1 In particular the 
^ • eigenvalues of the Laplace operator on a graph are typically distributed as 

the eigenvalues of large random matrices where the appropriate random 
matrix ensemble for comparison is determined by the symmetries of the 
quantum system; this correspondence is known as the Bohigas-Giannoni- 
Schmit conjecture. 2 

Quantum graphs were introduced as a model of quantum chaos by Kot- 
tos and Smilansky. 3 Graph models are also of current interest in many areas 
of mesoscopic physics like Anderson localization, photonic crystals, micro- 
electronics, nanotechnology and the theory of wave guides; see Ref. 4 for a 
review. In this article we present results for the vacuum energy, 5,6 spectral 
determinant 7-9 and heat kernel asymptotics 10,11 of graph Laplacians. The 
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unified approach we adopt, based on the spectral zeta function of the graph, 
provides new results in each case. The full description of the construction 
of zeta functions of quantum graphs will appear in Ref. 12. 

2. Graph Laplacian 

Let G = (V, B) be a graph where V is the set of vertices and B is the set of 
bonds (or edges); for an example see Fig. 1. Each bond b € B joins a pair of 
vertices and for notational convenience we assume that bonds are directed 
with an initial and terminal vertex specified via the functions o : B — > V 
and t : B — > V respectively. We insist that the set B is symmetric so that 
b e B if and only if there is another bond b € B the reversal of b such that 
o(b) — t(b) and t(b) = o(b). The total number of bonds B = \B\ is therefore 
twice the physical number of connections in the network. Each bond b is 
associated with an interval [0, Lb] so that the local coordinate Xb = at o(b) 
and Xb = Lb at t(b). We naturally require that L^ — Lb so there is a single 
length associated with a physical connection in the network, consequently 
x-fj = Lb — Xb- The total length of the graph G is denoted by C — | J2b=i Lb- 




Fig. 1. A star (or hydra) graph. 

A function ip on G is defined by the collection of functions {ipb{xb)}beB 
on the bonds of the graph. We require ^b{xb) = ^{Lb — £&) in order that 
ip be physically meaningful. The Hilbert space on the graph is 

L 2 (G)=0i 2 ([O,i b ]) . (1) 

6=1 

,2 

To define a Laplacian on G we consider the operator —-^i on the bonds 
with matching conditions between functions ^6 meeting at vertices of G to 
make the graph Laplacian self- adjoint: for example Neumann like matching 
conditions require ip to be continuous at the vertices and J2b\o(b)=v V'b(^) = 
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for all vertices v € V. All such self-adjoint matching conditions were 
classified by Kostrykin and Schrader 13 and we follow their scheme. The 
matching conditions on the whole graph are specified by the pair of B x B 
matrices A, B via 

Aip + Mip' = , (2) 

where ip = (V>i(0), . . . , ^b(0)) T and ip' = (tp[(0), . . . , ip' B (0)) T . These 
matching conditions define the domain of a self-adjoint Laplace operator iff 
rank(A,B) = B and AB 1 * = BAt. 

To find the graph spectrum we want to solve the Laplace eq. on the 
bonds 

d 2 

-jl^t 1 !-) = k 2 i> h (x b ) . (3) 
ax b 

In the case of a star graph, Fig. 1, solving Eq. (3) on the bonds and substi- 
tuting in the matching condition at the center one obtains the well known 
secular eq. of the star graph X^i tan kL b = 0, whose solutions {kj\kj ^ 0} 
are square roots of the eigenvalues of the Laplacian on the star. 

If instead one substitutes solutions of (3) in the general form of matching 
conditions (2) we find the sequence {kj} are the solutions of a secular eq. 
f(k) = where, 

/(*) = det ( A + m ( ~ C °^ L) CSC( ^ T } y\ , (4) 
3 w V V csc(fcL) - cot(fcL) J J ' w 

and L = diag{Li, . . . , Lb}- This can be compared to the scattering matrix 
formulation introduced by Kottos and Smilansky. 3,14 

3. Results 

The zeta function of the quantum graph can be formulated as a contour 
integral applying the argument principle to (4) following the technique in- 
troduced in Refs. 15,16, namely 

CM - 1 <*7* - 1 - h j^hw* ■ <» 

The prime on the sum denotes the exclusion of zero modes if present. Gener- 
ically the poles of / are the whole of the set {mit/L^m e Z, b = 1, . . . , B} 
and this is the case we consider. a The contour c is then chosen so as to 



a It is possible for particular choices of matching conditions A, B and bond lengths {Lf,} 
that individual poles cancel in the determinant, however, a small perturbation of the 
bond lengths removes this. 
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include the zeros of / on the positive real axis and exclude the poles. Mak- 
ing a contour transformation to an integral along the imaginary axis closed 
with a semicircular arc on the right we obtain an expression for £(s) valid 
in the strip < Sfts < 1, 



«s)=Ch(2s)^(-) + —J o t---log/(t)dt, 



(6) 



f{t) = det ( A - tB ( COth ^ -csch(tL) \ \ 

• /w V V-csch(tL) coth(tL) J J w 

where (r is the Riemann zeta function. 

Regularized formulations of the vacuum energy, and the spectral deter- 
minant as well as the heat kernel asymptotics are obtained directly from the 
zeta function. The vacuum energy, formally E = | J2jLo 'kj where fc| is an 
eigenvalue of the Laplacian, is defined as E c = ^C( — 1/2)- Differentiating 
with respect to the Lp, the length of the bond (3, we obtain the Casimir 
force on the bond, 



F? = t£t* + - —\o g f{t)At. (8) 



i r x d 

The spectral determinant is formally the product of the eigenvalues 
n^Lo 'kj- ^he zeta function representation (6) allows a direct evaluation 
of the regularized spectral determinant, as 

det'(-A) = cxp ( - C'(0)) = J B { ' . (9) 

Here, cat denotes the first non-zero coefficient in the t — > oo asymptotic 
expansion of 

f(t) ~ dot (A - tB) = det B t 2B + c 2B -it 2B ~ 1 H h C\t + det A . (10) 

So cat = detB when detB ^ 0. The factor of cat in (9) is introduced by the 
analytic continuation of ((s) to s = 0. 

The heat kernel for a quasi one-dimensional system has an expan- 
sion for t -> of the form, K{t) = e"^' ~ ££=0,1/2,1,... M £ ~ 1/2 - 
The heat kernel coefficients are related to the zeta function 17 by at — 
Res(C(s)r(s))| s= i/ 2 _£. Subtracting sufficiently many terms in the t — > oo 
expansion of /, the zeta function (6) can be extended as far left of the 
imaginary axis as required and we obtain the full asymptotic expansion of 
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the heat kernel of the general graph Laplacian, 

The coefficients b n are obtained from the expansion of log f(t) at infinity, 

OO 7 

log f(t) ~ N log t + log c N + £ -£ . (12) 

n=l 

For a given set of matching conditions defined by matrices A, B the coeffi- 
cients &„ can be easily obtained with a computer algebra package. 
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